Stochastic processes and conformal invariance 



(N 
O 

o 

(N 

> 
O 

(N 



o 



Jan de Gier 1 , Bernard Nienhuis 2 , Paul A. Pearce 1 and Vladimir Rittenberg 1,3 
1 Department of Mathematics and Statistics, University of Melbourne, Parkville, Victoria 3010, Australia 
2 Universiteit van Amsterdam, Valckenierstraat 65, 1018 XE Amsterdam, The Netherlands 
3 Physikalisches Institut, Bonn University, 53115 Bonn, Germany. 
(Dated: February 1, 2008) 

We discuss a one-dimensional model of a fluctuating interface with a dynamic exponent z = 1. 
The events that occur are adsorption, which is local, and desorption which is nonlocal and may 
take place over regions of the order of the system size. In the thermodynamic limit, the time 
dependence of the system is given by characters of the c = logarithmic conformal field theory of 
percolation. This implies in a rigorous way a connection between logarithmic conformal field theory 
and stochastic processes. The finite-size scaling behavior of the average height, interface width 
and other observables are obtained. The avalanches produced during desorption are analyzed and 
we show that the probability distribution of the avalanche sizes obeys finite-size scaling with new 
critical exponents. 
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The structure of growing interfaces continues to be a 
subject of major interest and a characterization of the 
various universality classes of critical behavior remains 
an open question [jlj. In this paper we present a one- 
dimensional adsorption-desorption model of a fluctuating 
interface which belongs to a new universality class where 
the dynamical critical exponent z = 1. In this model the 
interface evolves following nonlocal Markovian dynamics 
(for an adsorption-desorption model with local dynamics 
and z > 2 see || ) . The relaxation rules are such that one 
has avalanches with a long tail in their probability distri- 
bution function (PDF). The model belongs therefore to 
the self-organized criticality class (SOC) || ||]. What 
makes the model special is that the correlation lengths 
in the time and space directions are both proportional 
to the size of the system (this is not the case for other 
SOC models |, §). The single scale that exists in the 
system is therefore its size. Moreover, our model has the 
big advantage of being solvable. 

The Hamiltonian which gives the time evolution of our 
model is integrable. The finite-size scaling (FSS) limit of 
its spectrum can be obtained using the Bethe Ansatz j^j] 
and is given by characters of a c = Virasoro algebra. A 
logarithmic conformal field theory (LCFT) with c = i 
appears also in other domains of physics such as systems 
with quenched disorder and the quantum Hall effect ||, 
and possibly string theory p0| . Moreover pl| , the PDF 
describing the stationary state is related to combinato- 
rial aspects of the ice model defined on a rectangle with 
special boundary conditions. This is probably the rea- 
son why we are able to conjecture exact expressions for 
various quantities connected to the model even for finite 
lattices. 

We consider an interface on a one-dimensional lattice 
of size L + 1 (L — 2n even). The non- negative heights 
hi (i = 0,1,.., L) that specify the interface, obey re- 
stricted solid-on-solid (RSOS) rules, 



Alternatively, we can describe the interface using slope 
variables Sj = (/ij+i — hi-\)/2, (sq = sl = 0). The 
interface evolves through adsorption and desorption ac- 
cording to the following rules. Adsorption, which locally 
changes the height hi to hi + 2, takes place with a rate 
equal to 1 at a local minimum of the interface (s» = 0, 
hi < (hi— i, /ij+i)). The non-vanishing rates for desorp- 
tion can be understood using the notion of active seg- 
ments. A segment of a configuration with the end-points 
a and b is defined by the conditions: h a = hb = h and 
hi > h for a < i < b. We call a segment active if at least 
one of the boundary slopes s a or Sb is nonzero. In the 
desorption event, all the b — a — 1 heights hi contained 
in a segment decrease by two units (hi i— > hi — 2) with a 
rate 



5(s a -l) + 5(s b + l), 



(2) 



H+l 



±1, 



ho 



0, hi > 0. 
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where 5 is the discrete Kronecker symbol, the other 
heights remaining unchanged. In order to find which 
desorption events can take place, for each of the C' n — 
(2n)\/((n+ l)(nl) 2 ) configurations one first looks at how 
many active segments one has and then one uses (^|). If 
the two slopes s a and Sb are zero, the desorption rate 
is zero. This observation has two consequences. First, 
in the stationary state we expect to see predominantly 
configurations with large terraces, i.e. intervals where 
the slope is zero for all the sites. Next, it is meaningful 
to consider clusters. Those are segments where the end- 
points have the heights equal to zero. According to the 
rules (|J) desorption takes place only within a cluster. 

The dynamics of the interface can easily be visualized 
using tiles (tilted squares) which cover the area between 
the interface and the substrate (hn = 0, h-a+i = 1)- 
For a given configuration c, the number of tiles is u(c) = 

\ Y^t=2 hi +1 — L/4:. According to the rules given above, 
through adsorption one adds one tile, through desorption 
one loses a layer with b — a — 1 tiles (this is an odd num- 
ber) . If one looks at our model from the point of view of 
self-organized critical phenomena, the number of tiles u 



2 



lost in a desorption event defines the size of an avalanche. 

The physics of our model can be understood in the 
following way. We have a gas of tiles in the presence of a 
surface. Tiles uniformly try to attach to the interface and 
are succesful on sites that are lower than their neighbors. 
Tiles get reflected back into the gas on sites that are 
higher than their neighbors and trigger desorption events 
and are reflected on local slopes. Because the stationary 
state mainly consists of terraces, desorption events occur 
less frequently, but when they occur, they can take larger 
number of tiles. 

Before describing the physical properties of the model 
we will give another formulation which enables us to ob- 
tain the FSS spectrum of the Hamiltonian, defined by the 
rates given above and thus to find the dynamic exponent. 
To each RSOS configuration with L+l sites one can as- 
sociate a configuration of n non- intersecting half- loops on 
L = 2n sites. The height hi is the number of half-loops 
above the midpoint of the sites i and i + 1 in the half- 
loop configuration. For example, the following picture 
illustrates the association between a configuration with 
three half-loops and an L = 6 RSOS configuration with 
two clusters and one tile. 



(3) 
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The half-loop configurations on the other hand, can 
be associated with the left ideal of the Temperley Lieb 
(TL) algebra at Q = 1 |l2| . The L — 1 generators of the 
algebra satisfy the relations 

ej = ei, eiei±iet = e i7 [e h ej] =0, \i - j\ > 1, (4) 

and have the graphical representation 



e, 



(5) 



-1 i i+1 i+2 



L-l L 



The left ideal is generated by the action of TL generators 
on Iq = Yli—i e 2i-i- For example the action of on the 
half-loop configuration (||) is 



(6) 



The Hamiltonian gives the time evolution in the vector 
space of the half-loop (RSOS) configurations, 




H 



L-l 

Ed 



Y^HMt). (7) 



P c (t) is the (unnormalized) probability to find the system 
in the configuration c at the time t and r c( i = —H c( i give 
the rates of our model for the transitions d — > c. Since 
H is an intensity matrix (Y] c ^ cd = 0) PI' ^ nas a zero 



eigenvalue with a trivial bra and a nontrivial ket which 
gives the probabilities in the stationary state, 



(0|F = 0, 
H\0) = 0, 



I = (1,1,- 



P c = lim P c (t). 

t—>oo 



(8) 



Let us compute the dynamic critical exponent z. We 
use the following representation of the TL algebra 



1 1 

4 ~ 2 
+ h(a, 



u i+i 

T i+1 



(9) 



where a x ,a v , a z are Pauli matrices, — 2A = q + q^ 1 = 1, 
2h = q — q^ 1 = iy^3 and q = e 7 "/ 3 . In this representation 
H given by (0) becomes the Hamiltonian of the XXZ 
quantum chain with L sites and the energy gaps Ek 
(the ground-state energy is zero for any L) scale like 



E k = nv s x k L z , 



(10) 



with z = 1, v s = 3\/3/2 is the sound velocity and 
Xfe are related to surface exponents of a c = LCFT of 
percolation In the continuum limit, the spectrum 

of the Hamiltonian (]7]) , and thus the values of x n , in the 
full TL algebra is given by the Virasoro characters [jl4j 



(11) 



fc=i 



Here w = exp(irTv s / L) parametrizes the temporal (T) 
and spatial (L) extent of the stochastic process. In the 
subspace of the RSOS configurations the spectrum is 
given by Xo(^)- To our knowledge, for the first time, 
a connection is made between stochastic processes and 
LCFT. This implies among other things that the space 
and time correlation lengths are the same and that, in the 
continuum limit, the forms of the space-time correlation 
functions in the stationary state are known |15|. 

As discussed in detail in 1 1 1 the stationary state of 



our model with L sites is related to the two-dimensional 
ice model |l6| defined on a rectangle of dimension L x 
(L — 1)/2 with special boundary conditions |l7|, [l8| . This 
model is equivalent to a fully packed loop (FPL) model 
p9[ on the rectangle, all configurations being equally 
probable. We briefly explain this connection and its con- 
sequences. If we choose the (unnormalized) probability 
of the "pyramid" configuration [s\ = ... = s n -i = L 
s n = 0, L = 2n) to be equal to one (this configuration 
has the smallest probability) , then the normalization fac- 
tor of the stationary state p0|] is equal to 



(0|Q) 



J=0 



(2j + l)!(6j + 3)! 
(4j + 2)!(4j+3)! : 



(12) 



This is precisely Ay(2n +1), the number of configura- 
tions of the FPL model on the rectangle. Moreover, the 
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(unnormalized) probabilities for the other half-loop con- 
figurations are integer numbers which are equal to the 
numbers of configurations of the FPL model with the 
external occupied edges connected in the same way as 
in the configuration of half- loops 111], E]J. This implies 
that the PDF describing the stationary distribution in 
the RSOS configurations space corresponds to a uniform 
probability distribution in the space of the FPL model, 
or equivalently the ice model. The RSOS configuration 
with the largest (unnormalized) probability is the one 
with the longest terrace (s± = .. = S£ = 0) in which the 
interface coincides with the substrate. The normalized 
PDF is obviously p c = P c /(0|0). 

We now describe some properties of the model. They 
were obtained either numerically for lattice sizes up to 
L = 18 or an exact expression was conjectured and 
checked up to L = 18. We start with the properties of 
the stationary state. The fraction of the interface covered 
by terraces [J22| 



is 



r(L) = j^£{l-\ S] 



3L 2 -2L + 2 
(L- l)(4L + 2)' 



(13) 



This implies that for large L, three quarters of the sur- 
face is covered by terraces. We now consider the number 
of clusters per configuration C — Y^j^i^i^j)- av- 
erage and the average of its square have the following 
expressions 



(C) 



1 ^ (2j 



n 



l)(3j + 4) 



3 to (J + I)(6j + 1) 



(C* 2 ) w (1.29 - 0.49L~ a63 )(C) 5 



0.738L 2/3 , (14) 
(15) 



From ( |l5|) which we deduce that the compressibility Q 
of the gas of clusters k ss 0.29L diverges like the size of 
the system. 

The average height (h) and interface width Woo, which 
characterizes the roughness of the surface in the station- 
ary state, 



1 L 



w 



(h 2 



(16) 



are compatible with the following behavior, 

(h) ss 0.141n(L/2) , Woc S3 0.34 (lnL/2) 1/3 . (17) 

This implies that the PDF of the heights has a very small 
dispersion and also that the usual exponents Q a = /3 
vanish. As in any fit involving logs, the formulas given in 
( |l7| ) are probably not the last word. What is certain is 
that the width grows slowly with the size of the system 
implying that the surface is only marginally rough. It 
is interesting to mention that marginally rough surfaces 
(with z = 1.581 corresponding to the directed percolation 
universality class) were also encountered at a critical 
point dividing a moving rough KPZ phase from a smooth, 



massive phase. In these models a factor 3 between the 
powers of the logs (like in (|l7j)) was also seen |25|. The 
factor 3 is typical of a large class of growth problems Q . 
To complete the interface growth picture, we calculated 
the ratio of the time dependent interface width w(t, L) 
to Woo for various lattice sizes. The curves were obtained 
by solving on a computer the differential equations (j?]) 
taking the substrate as the initial configuration. If z = 1 
and if the Family- Vicsek scaling applies, this ratio should 
converge to a function oit/L. We have checked that this 
is indeed the case. 

We now consider the avalanches in our model. If in the 
stationary state one considers the configuration d with 
u(d) tiles, it changes with rate r c d into the configuration 
c with u{c) tiles. The rate of changes in which u tiles are 
lost, can be written as 



R(u, L) = 2_, 3(u(d) — u(c) — u)r c 

c^td 



(18) 



We have found the following values for the rate of ad- 
sorption of one tile, R(u = — 1,L), and the rate of all 
possible avalanches, R(u > 0,L), 



R(-1,L) 



3L(L-2) 



(19) 



4(2L + 1) 

R(u>0,L) = L-2-2R(-l,L). (20) 

Notice that for large L the rate for the adsorption of one 
tile (u = —1) is 3/2 times larger than the rate for an 
avalanche (u > 0). This explains the relative rarity of 
the desorption events. 

Since through desorption one loses an odd number of 
tiles, it is convenient to write u = 2v — 1 and to consider 
v as the size of an avalanche. Given the occurrence of an 
avalanche, its size v is distributed according to the PDF 



S(v,L) = 



R(2v-1,L) 



(21) 



R{u > 0,L) ' 

FSS U predicts the following form for this PDF, 

S{v 7 L)=v- T F{v/L D ). (22) 

One way to get the exponents of the FSS function is to 
consider the moments IM, 



{v k ) =]r« fe S(«,L)~L 



a(k) 



v=l 



and one expects 



a(k) 



0, 

D(k + 1 -t) 



k < t 
k > t 



(23) 



(24) 



From equation (|20|) one gets (v) = (5L + 4)/4(L + 2), 
which implies that for large L the average size of an 
avalanche is 3/2 tiles. A numerical investigation of the 
other moments for 1 < k < 5 indicates cr(k) — k — 2. 
This implies D = 1 and r — 3. The numerics can be 
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FIG. 1: Avalanche scaling function F(v/L). The data are 
obtained for v > 1 and L = 8, .... 18. 



quite precise even if we have data for L up to 18. This 
is because we know S(v, L) exactly and we can use VBS 
approximants [^8| to derive the large L behavior of the 
moments (the convergence is less good around tr(fc) = 
due to logarithmic effects). The value D = 1 was to be 
expected since L is the only characteristic length in our 
system. A consistency check was done assuming D = 1 
in (^2|) to see for which value of tau one gets data collapse 
for the scaling function F(v/L). Since we have data up 
to L = 18 only we cannot expect a precise value of r nor 
of F. Nevertheless as shown in Fig. H, a data collapse is 
visible for r = 3.2. 



To conclude, we have presented an SOC model of a 
critical fluctuating interface belonging to a new univer- 
sality class. We have also shown the connection between 
the model and a LCFT with c = 0, which among other 
things implies that the dynamic exponent z = 1. Several 
FSS exponents of expectation values have been computed 
in the stationary state but their identification with the 
scaling dimensions of the LCFT has not been completed. 
The FSS exponents of the avalanche PDF have also been 
determined. Much is still to be understood in this model, 
such as correlation functions and its off-critical behavior. 
We hope to come back to these topics in a future publi- 
cation. 
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